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ABSTRACT 



Recently, the WMAP group has published their five-year data and consid- 
ered the constraints on the time evolving equation of state of dark energy for the 
first time from the WMAP distance information. In this paper, we study the 
effectiveness of the usage of these distance information and find that these com- 
pressed CMB information can give similar constraints on dark energy parameters 
compared with the full CMB power spectrum if dark energy perturbations are 
included, however, once incorrectly neglecting the dark energy perturbations, the 
difference of the results are sizable. 

Subject headings: Cosmology: theory — (Cosmology:) cosmic microwave back- 
ground — (Cosmology:) cosmological parameters 



1. Introduction 



The newly releas e d Wilkinson Micr o wave Anisotropy Probe five-year data (WM AP5) 
jHinshaw et al. 20081 : iNolta et al. 20081 : bunklev et al. 20081 : iKomatsu et al. 2008ah . de- 
tecting the Cosmic Microwave Background (CMB) to an unprecedented precision, make 
it possible to improve the constraints on almost all the cosmological parameters, including 
the equation-of-state (EoS) w of the unknown energy budget, dark energy. Defined as the ra- 
tio of pressure over energy de nsity, w = P/p, EoS can be used to clas s ify various dark energy 



models, such as quintessenc e (jWetterich 19951 : iRatra & Peebl es 1998: JPeebles fc Ratra 1998 
Wetterich lQgi), phantom JCaldwell 2002|), quintom (IFeng et al. 20051) k-essence 



( Armendariz-Picon et al. 200ol : Armendariz-Picon et al. 2001 ). etc., which is of great the- 
oretical significance to unveil the mystery of dark energy. Therefore, trying to study the 
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evolu tion history of EoS of dark energy plays a crucial ro l e in modern observ ational cosmol- 
ogy (IHuterer fc Starkman 20031 ; IWang &: Tegmark 20051 ; IZhao et al. 2007a ). Simply put, 
one can choose an arbitrary parametrization of w and constrain the introduced dark energy 
parameters from the astronomical observational d ata, including CMB, Supernova type la 
(SN la), Large Scale Structure (LSS) and so forth JWang fc Mukheriee 20071 ; IWright 20071 : 



Zhao et al. 2007al : I Zhao et al. 2007bf ). 



Recently, the WMAP group has released their five-year data and for the first time 
considered the constraints on the time evolving EoS of dark energy using the WMAP distance 
information. This method has the advantage of reducing computation time by orders of 
magnitude, yet the effectiveness and the level of approximation compared with the full CMB 
power spectrum computation remain unclear. In this paper, we make a thorough test of 
this simplified method to investigate whether it is safe to constrain dark energy with time 
evolving EoS. Our paper is structured as follows: In Section II we describe the method and 
the data; in Section III we present our main results; finally we present our conclusions in 
Section IV. 



2. Method and data 



To study the dynamical behavior of dark e nergy, we choose the parametrization of 



the time evolving EoS of dark energy given by (jChevallier fc Polarski 200ll : iLinder 20031 : 
Komatsu et al. 2008al ): 

w(a) =w + w a (l - a) , (1) 

where a = 1/(1 + z) is the scale factor and w a = —dw/da characterizes the "running" of the 
EoS (RunW henceforth). For the ACDM model, Wq = — 1 and w a = 0. 

When using the MCMC global fitting strategy to constrain cosmological parameters, 
it is crucial to include dark energy perturbations, especially for the time evolving EoS 
of dark energy models. This issue h as been realized by many researchers including the 



WMAP group dWeller fc Lewis 20031 : lYeche et al. 2009 : IZhao et al. 20051 : Kia et al. 2006 



Spergel et al. 20071 ). However, when the parameterized EoS crosses —1, one cannot han- 
dle the dark energy perturbations based on quintessence, phantom, k-essence and other 
non-crossing models. By virtue of quintom, the perturbations at the crossing points are con- 
tinuous, thus we have proposed a technique to treat dark energy perturbations in the whole 
parameter space. For details of t his method, we refer the readers to our previous papers 
Jzhao et al. 20051 : kia et al. 2006h . 



In this study, we have modified the publicly available Markov Chain Monte Carlo pack- 
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age CosmoMC0 ( Lewis &: Bridle 2002 ) to include the dark energy perturbations with EoS 



across —1. Furthermore, we assume purely adiabatic initial conditions and a flat universe. 
Our most general parameter space is: 

P = (u b ,u c , Q s ,r,w ,w a ,n s ,\n(10 10 A s )) , (2) 

where io b = VL b h 2 and uo c = Q c h 2 , where fit, and Q c are the baryon and cold dark matter 
densities relative to the critical density, Q s is the ratio (multiplied by 100) of the sound 
horizon at decoupling to the angular diameter distance to the last scattering surface, r is 
the optical depth to reionization, A s and n s are the amplitude and the tilt of the power 
spectrum of primordial scalar perturbations. For the pivot scale of the primordial spectrum 
we set K = 0.05 Mpc" 1 . 

The WMAP distance information used by WMAP group include the "shift parameter" 
R, the "acoustic scale" I a and the photon decoupling epoch z*. R and I a correspond to 
the ratio of angular diameter distance to the decoupling era over Hubble horizon and sound 
horizon at decoupling respectively, given by 



R(z*) = y/n m H* X (z.) , (3) 
l A (z m ) = kx{z*)IXs{z*) , (4) 

where x( z *) an d Xs{z*) denote the comoving distanc e to z* and comoving s ound horizon at 



z* respectively. The decoupling epoch z* is given by (jEisenstein &: Hu 19981 ) 

z, = 1048[1 + 0.00124{Q b h 2 )-°- 738 }[l + gi (fl m h 2 ) 92 } , (5) 

where 

_ 0.0783(f^ 2 )-°- 238 _ 0.560 
91 ~ l + 39.5(fi 6 /i 2 ) - 763 ' 92 ~ l + 21.1(n 6 /i 2 ) 1 - 81 ' 

The WMAP distance information encode in part of the CMB information and can 
constrain cosmo logical parameters to some extent. It is worth carefully investigating on the 
effectiveness of the constraints from the distance information compared with the full CMB 
power spectrum computation. To do this, we follow the procedure shown in the flow chart: 



Full WMAP5 Data M ™ C w ,w a 

MCMC | $ Compare 

7 D MCMC /„s 

l A ,R,z* -> w ,w a (7) 



http : / / cosmologist . info /cosmomc / . 
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which are detailed as follows: 



1. Making a global fitting with MCMC method to constrain wq, w a , and also I a, R and 
z* using the full WMAP5 power spectrum. In this step, we have done two types of 
calculations, one with and the other without dark energy perturbations; 

2. Using the resultant I a, R and z* to constrain dark energy parameters wo and w a ; 

3. Comparing the results of the constraints on wq and w a obtained from step 1 and 2. 



In step 1, we calculate the likelihood of CMB power spectrum using the routine supplied 
by the W MAP grourH In step 2, we calculate the likelihood of WMAP distance information 



as follows (IKomatsu et al. 2008al ): 



X 2 = -2\nL = {xf - xf ata )(C- 1 \ 



i th _ „data\ 



where x = (Ia,R,z*) is the parameter vector and (C is the inverse covariance matrix 
for the WMAP distance information. 

Since the purpose of this paper is not to make a global analysis, in order to see the eff ects 
of the other cosmological data, we include the gold sample of 182 SN la (IRiess et al. 20071 ) for 
a joint constraint on EoS of dark energy in the combination with the WMAP5 data. In this 
study we also make use of Hubble Space Telescope (HST) measurement of the Hubble param- 
eter H = 100/i km s _1 Mpc -1 by multiplying the likelihood by a Gau ssian likelihood functio n 
centered around h = 0.72 and with a standard deviation o = 0.08 fjFreedman et al. 20011 ) . 
and a Gaussian prior o n the baryon dens ity VL^h 2 = 0.022 ± 0.002 (la) from Big Bang 
Nucleosynthesis (BBN) jBurles et al. 200lh . 



3. Results 

The WMAP distance information are extracted from the full WMAP5 power spectrum 
by assuming a certain cosmological model, and they should be model dependent. In Fig JT] we 
present the one dimensional distributions of the WMAP distance information for different 
cosmological models. 

In the upper three panels of Fig{T] we show the distributions of I a, R and z* for five 
cosmological models: flat ACDM model; ACDM with curvature; flat ACDM model with 
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massive neutrinos, with running of spectral index and with tensor perturbations, respectively. 
We find that the distributions of R and are quite different in these five cases, while the 
acoustic scale I a does not change significantly. These results indicate that when using these 
distance information to constrain cosmological parameters, one should be clear about the 
assumed cosmological model. In Table I we also list the median la constraints on the WMAP 
distance information from the full WMAP5 data for different cosmological models. 




300 302 304 1.65 1.7 1.75 1088 1090 1092 1094 



Fig. 1. — One dimensional posterior distributions of Ia, R and z* with the WMAP5 data 
for different cosmological models. In the upper panels, the black solid line is given by the 
standard flat ACDM model, while the red dashed line, the blue dash-dotted line, purple 
solid line and the green dotted line are given by ACDM with non-zero fl*., flat ACDM with 
f v , with a s and with r respectively. In the lower panels, the black solid line is still from the 
standard ACDM model, while the red dashed line and the blue dash-dotted lines are given 
by the dark energy model with time evolving EoS (RunW model) with and (incorrectly) 
without dark energy perturbations respectively. 

In the lower three panels of Fig{T]we show the results for three flat models with different 
dark energy properties: ACDM model, RunW model with and (incorrectly) without dark 
energy perturbations. These results do not show significant differences in the WMAP dis- 
tance information among different dark energy models. We also compare the results obtained 
with and (incorrectly) without dark energy perturbations and find that simply switching off 
dark energy perturbations does not bias the results much at this stag^f). In the following 



3 The distance information are determined by the background parameters and not affected by dark energy 
perturbations significantly. 
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calculations, we use the WMAP distance information obtained from the RunW model with 
dark energy perturbations included. The corresponding inverse covariance matrix is shown 
in Table II. 




Fig. 2. — 68% and 95% confidence levels constraints on (wo,w a ) from full WMAP5 data 
and WMAP distance information respectively. Red solid lines are obtained from the full 
WMAP5 data including dark energy perturbations; black dash-dotted lines are from the full 
WMAP5 data incorrectly neglecting dark energy perturbations; and blue dashed lines are 
from WMAP distance information. 

In Figf2]we compare the constraints on wq and w a obtained from the full WMAP5 power 
spectrum with the one obtained from WMAP distance information given in Table I and II. 
From this plot we can see that the WMAP distance information and the full WMAP5 power 
spectrum with dark energy perturbations included can give quite similar constraints on wq 

Table I. Median ler constraints on WMAP distance information using full WMAP5 data for different cos- 
mological models. 



Models 


U 


R 


z* 


ACDM 


302.15 ±0.842 


1.71 ±0.021 


1090.92 ±0.969 


ACDM + n K 


302.32 ±0.899 


1.72 ±0.021 


1091.26 ± 1.004 


ACDM + m v 


302.30 ±0.873 


1.75 ±0.031 


1091.98 ± 1.244 


ACDM + a s 


302.36 ±0.878 


1.74 ±0.031 


1092.72 ± 1.817 


ACDM + r 


301.76 ±0.944 


1.69 ±0.027 


1089.72 ± 1.366 


RunW with pert. 


302.20 ±0.865 


1.72 ±0.021 


1091.10 ±0.991 


RunW w/o pert. 


302.14 ±0.875 


1.71 ±0.021 


1090.97 ±0.985 
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and w a . However, when the dark energy perturbations are incorrectly switched off (black 
dash-dotted lines in Fig|2j), the results between the two methods are quite different. 

The WMAP distance information mainly include the information of the oscillatory struc- 
tures of the CMB power spectrum, which come from the small angular scale (large /) of the 
power spectrum. On the other hand, for the full CMB power spectrum, they combine more 
information than the distance information, especially at large angular scale (small I) . At large 
angular scale, they are affected by the late Integrated Sachs- Wolfe (ISW) effects, which are 
dark energy dependent. Thus tighter constraints on (w , w a ) are anticipated by using the full 
CMB spectrum than those from using the distance priors only. This is clearly demonstrated 
in Figf5] (dashed contours versus dash dotted one). It is noted that the dash dotted contours 
are calculated without including dark energy perturbations, and thus the constraining power 
of the late ISW effect on dark energy parameters is fully realized. However, when including 
the dark energy perturbations, which are mainly effective at small I, the constraints on dark 
energy parameters from the late ISW effects are signific antly reduced, re sulting in similar 



contours shown by the dashed and solid lines in FigOu ( jKomatsu 2008bl ). The differences 
between the solid and the dash dotted contours also show that how biased results can be 
obtained if the dark energy perturbations a re incorrectly neglected in the full CMB data 



analysis (jXia et al. 20061 ; ISpergel et al. 20071 ) 



In Figj3j we give the constraints on dark energy parameters by adding the SN la data. 
We can see that the constraints on dark energy parameters are tightened and the differ- 
ences between the results obtained from "full WMAP5 power spectrum + SN la" and from 
'WMAP distance information + SN la" become insignificant. 



In our analysis we use a specific parametrization the "RunW" model. We expect that our results 
hold qualitatively for other dark energy parameterizations. Quantitatively however, the specific results are 
dependent on the detailed calculations on different dark energy parameterizations that are used. 

Table II. Inverse covariance matrix for the WMAP distance information I a, R and z* in 
RunW dark energy model when including dark energy perturbations. 







R(z*) 






1.795 


31.596 


-1.146 


R(z*) 




5409.68 


-94.58 








2.891 
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Fig. 3. — Two-dimensional constraints on the parameters of dark energy from the combined 
WMAP5 and SN la. The red solid and black dotted lines are the constraints from "full 
WMAP5 data + SN la" with dark energy perturbations and "WMAP distance information 
+ SN la" respectively. 

4. Summary 

In this paper, we have studied the effectiveness of the WMAP distance information on 
constraining the dark energy parameters, by comparing with the full WMAP5 power spec- 
trum analysis. We first present the level of the model dependence of the distance information 
in different cosmological models. We further clarify that by taking into account dark energy 
perturbations properly, the WMAP distances can give unbiased information on dark energy 
parameters relative to the full CMB analysis. 
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